Observation of twin beam correlations and quadrature entanglement by frequency 

doubling in a two-port resonator 
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We demonstrate production of quantum correlated and entangled beams by second harmonic 
generation in a nonlinear resonator with two output ports. The output beams at A = 428.5 nm 
exhibit 0.9 dB of nonclassical intensity correlations and 0.3 dB of entanglement. 
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Continuous light beams that exhibit nonclassical 
statistics are of interest as a tool for studying quan- 
tum fieldsQ and for a number of applications that in- 
clude precision measurements 0, writing subwavelength 
spatial structures'^, and as resources for quantum infor- 
mation and communication protocols (Jj. The most suc- 
cessful and widely used approach to generating nonclassi- 
cal light employs parametric down conversion in crystals 
with a quadratic nonlinearity. At the microscopic level 
non-classical correlations and entanglement arise due to 
the possibility of converting a single high frequency pho- 
ton at 2uj into a pair of entangled lower frequency photons 
at uj. 

In order to experimentally generate nonclassical beams 
with carrier frequency uj one typically starts with a co- 
herent source at lu which is frequency doubled to 2w. The 
light at 2uj is then used to drive a downconversion pro- 
cess to generate nonclassical light at frequency lu. These 
multiple steps add to the complexity of the experimen- 
tal arrangement and limit the possibility of generating 
nonclassical light at high frequencies. In this letter we 
demonstrate for the first time that nonclassical intensity 
correlations, as well as quadrature entangled beams, can 
be generated directly by frequency upconversion. In this 
way a coherent source at frequency lu produces quantum 
correlated beams at frequency 2uj. 

Consider the interaction geometry shown in Fig. ^ 
where a beam of frequency lu pumps a resonator that 
has two exit ports for the second harmonic beams at 
frequency 2a;. The cavity mirrors are assumed perfectly 
transmitting for the harmonic beams which are generated 
in a single pass of the intracavity pump field through the 
nonlinear crystal. It is well known that second harmonic 
generation (SHG) results in squeezing of the fundamen- 
tal and harmonic beams |^. The generation of multi- 
beam correlations in second harmonic generation is less 
well studied than in the case of parametric down con- 
version. Calculations have demonstrated the existence 
of correlations between the fundamental and harmonic 
fields including entanglement between the fundamen- 
tal and harmonic fields |2j and entanglement in type II 
SHG in the fundamental fields alone[S|. The possibility 
of nonclassical spatial correlations in either the funda- 
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FIG. 1: (color online) Conceptual schematic for generation of 
intensity correlations and quadrature entanglement by SHG 
in a two-port resonator. 



mental or harmonic fields alone and of entanglement 
in the fundamental field, IQj has also been shown in mod- 
els that include diffraction. Here we consider a situation 
where there are two harmonic output beams that share 
the same intracavity pump field. The common pump 
field couples the output beams, and our recent analysis 
of this interaction geometry^3| reveals that the two out- 
puts exhibit nonclassical intensity correlations, as well as 
quadrature entanglement. 

Nonclassical intensity correlations can be demon- 
strated by comparing the fluctuations of the sum and 
difference of the photocurrents generated by the output 
beams. We combine the detector outputs to form the 
quantities Id/s = Ii T 9I2 where g is an electronic scale 
factor. Calculations that account for the propagation of 
fluctuations in the two-ported cavitvprj show that with 
optimal choice of the scale factor g the variance of the 
sum and difference photocurrent fluctuations is 
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Here Sx is the amplitude quadrature noise spectrum of 
each harmonic output, Cx is the correlation coefficient 
of the amplitude quadratures at the two output ports, 
and we have simplified the more general expression given 
in Ref. to the case where the two output beams 

have identical squeezing spectra. It can be shown that 
the difference photocurrent is extremely close to the shot 
noise level while the noise of the sum photocurrent is 
suppressed below the shot noise. The reduced noise in 
the sum photocurrent arises from the amplitude squeez- 
ing of each beam and from the quantum correlations 
of the beams generated by the two-port cavity. As we 
demonstrate experimentally below the noise reduction of 
the intensity sum is greater than that due to the ampli- 
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tude squeezing of each beam alone which proves the role 
played by the shared intracavity pump beam in creating 
a nonzero correlation Cx- 

Calculations show that the output beams are entangled 
according to the Duan et al. inseparability criteria|l2j 
and exhibit EPR correlations for sufficiently high pump 
power. It can be shown that the beams are inseparable 
when the inequality[i3 V = \{2Sx+'2Sy+Cx-Cy) < 1 
is satisfied where Sy is the phase quadrature noise spec- 
trum of each beam, and Cy is the phase quadrature 
correlation coefficient of the beams. We see that veri- 
fying the presence of quadrature entanglement between 
the beams is technically more difficult than intensity cor- 
relation measurements since it requires determining both 
the amplitude and phase fluctuations. As the beams gen- 
erated by SHG have a nonzero mean amplitude, mea- 
surement of the phase quadrature using a strong local 
oscillator is inconvenient due to detector saturation. In- 
stead, we use an entanglement witness demonstrated by 
Glockl et al.0, 0| which allows us to verify quadra- 
ture entanglement using only intensity measurements. 
The beams are combined on a 50/50 beamsplitter with a 
tt/2 phase difference and the intensity fluctuations of the 
output beams are measured. The variances of the sum 
and difference photocurrents after the beamsplitter are 
(A|i+|)2 = 2Sx + Cx and (A|i_|)2 = 2Sy - Cy . The 
inseparability criterion can then be expressed as 
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which can be verified on the basis of intensity measure- 
ments alone. 

We have observed nonclassical beam correlations in a 
dual-ported resonator as shown in Fig. [3 The experi- 
ment uses an Ar-ion pumped Ti:Sa laser operating at 857 
nm which is frequency doubled with a KNbOa crystal in 
a confocal linear resonator. In order to suppress low fre- 
quency intensity noise the beam from the Ti:Sa laser was 
locked to a mode-cleaning cavity with linewidth 1.5 MHz 
(FWHM) usinga standard radio frequency (rf ) modula- 
tion technique[l5j. The beam was then mode-matched 
to a dual ported linear cavity with two plano-concave 
mirrors with radius of curvature R = 10 mm contain- 
ing a 1 cm long a-cut KNbOa crystal with anti-reflection 
coated ends. The input mirror had Tgsynm = 4% and 
ri-428nm < 5 %. The output mirror had RssTnm > 99 % 
and R428nm < 5%. Hence, only the fundamental field 
was resonant in the cavity. The separation of the mir- 
rors was set to 15.6 mm for confocal operation and the 
low power cavity finesse was measured to be 120. Using 
the measured finesse and the input coupler transmission, 
the roundtrip intracavity loss was calculated to be about 
1.1 %. 

Phase matching was controlled by varying the temper- 
ature of the crystal to produce two counter-propagating 
bright blue beams of wavelength 428.5 nm. The tem- 
perature tuning curves for the two harmonic outputs are 




FIG. 2: (color online) Experimental setup for entanglement 
generation with a dual ported resonator. FI = Faraday iso- 
lator, EOM — electrooptic modulator, DM = dichroic mir- 
ror. By inserting the mirrors shown with dashed lines de- 
tectors Di , D2 are used to verify intensity correlations. De- 
tectors Da,Dt are used for entanglement measurements. All 
detectors are Hamamatsu S3590-19 with quantum efficiency 
of 7? = 0.96 at 429 nm. 



shown in Fig.|3| The second output (counterpropagating 
to the incident pump beam) resembles the standard SHG 
phase matching curve for propagation without a cavity 
while the first output (copropagating with the incident 
pump beam) shows two distinct maxima. We believe the 
shape of the curves and the difference in optimal phase 
matching temperatures is due to a small but finite reflec- 
tivity of the cavity mirrors at the harmonic wavelength 
which changes the intracavity boundary conditions. A 
related sensitive dependence of the relative phase of the 
fundamental and harmonic fields on cavity parameters 
was noted in our earlier studies of self-pulsing in SHG [3. 
The experiments described below were performed using a 
crystal temperature of about 16.8 C which equalizes the 
output powers as indicated in the figure. In this case we 
expect the theoretical model used in Ref. [iJl to provide 
an accurate description of the experiment. 

The SHG cavity was locked to be resonant with the 
laser frequency uj to generate stable harmonic outputs 
in both directions. At Ppump = 34 mW incident on the 
cavity, 8 mW of 428.5 nm light was generated in each 
output beam. By inserting movable mirrors in the output 
beams detectors Di and D2 were used to measure the 
intensities of the two outputs as shown in Fig|21 The 
detector outputs were amplified and then recorded with 
a dual channel analog to digital converter (CompuScope 
14200, sampling frequency = 200 Msamples/sec, 14 bit 
data) . The noise power spectra of the difference and sum 
of the photocurrents were computed using a standard 
Fourier transform technique. 
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FIG. 3: (color online) Measured phase matching curves for 
the first (solid line) and second (dashed line) outputs from 
the dual ported resonator at Ppump ~ 31 mW. The phase 
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change of about 0.25 C. 



TT corresponds to a temperature 




7.5 10 12.6 
Frequency {MHz) 

FIG. 4: (color online) Spectra of the intensity noise for the 
intensity difference and sum of the two second harmonic fields. 
The broad response peak centered at 3.5 MHz is due to the 
filtering circuit used for AC-coupling of the photocurrent. The 
sharp peak at 15.8 MHz is due to residual modulation from 
the rf cavity lock. The data was acquired in 80 msec and the 
effective numerical resolution bandwidth is 100 kHz. 



The measured noise spectra of the intensity sum and 
difference correlations with g — 1 are given in Figs.^jand 
|S| The shot noise limit (QNL) for the intensity correla- 
tion is given by the noise spectrum of the intensity differ- 
ence of the two output ports. The maximum non-classical 
intensity correlation (corrected for electronic noise) was 
-0.90 ± 0.15 dB at a noise frequency of 6 MHz. The 
experimentally measured and theoretical noise reduction 
calculated from Eq. agreed to better than 5 % over 
the range of about 4 — 8 MHz. At lower frequencies the 
effect was masked by technical laser noise and at higher 
frequencies roll-off in the detector response reduced the 
magnitude of the observed effect. The amplitude noise 
of one of the output beams is shown by the blue line. We 
see that the noise reduction in the intensity sum clearly 
exceeds the amplitude squeezing of the individual beams 
which establishes the presence of the additional correla- 
tion Cx between the beams created in the SHG cavity. 
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FIG. 5: (color online) Noise spectra of the intensity sum(red 
curve) and amplitude squeezing of the first harmonic output 
measured at an output power of 9 mW(blue curve). The am- 
plitude squeezing of the second harmonic output agrees with 
the first output to ±0.1 dB. The shot noise reference level 
(QNL) was determined for the amplitude squeezing measure- 
ments using the difference of two detectors placed in one out- 
put port. The theory curves for the intensity sum and squeez- 
ing(dashed lines) were calculated using the cavity parameters 
given in the text and with a single pass conversion efficiency 
of 0.0059 W"^ which was determined from the mea- 

sured harmonic output power. The experimental data were 
corrected for electronic noise, but no correction was made for 
detector quantum efficiency and optical losses between cavity 
and detectors estimated to be 5 %. 



In order to verify if the two-beam correlations are 
strong enough to provide entanglement we measured the 
sum and difference intensities after mixing on the 50/50 
beamsplitter shown in Fig. |21 The first mirror on the 
right side of the SHG cavity was mounted on a piezo to 
allow active correction of the path length difference be- 
tween the two interfering arms. The photocurrents from 
Da and D;, were recorded with a dual channel analog to 
digital converter and analyzed in the same way as for the 
measurements of the intensity correlations. In this con- 
figuration the difference of the photocurrents from Di 
and D2 gives the corresponding shot noise level of the 
two entangled beams The shot noise level was cal- 
culated as (A|ii2|)^ = (A|zi — 0.95z2|)^ where the gain 
factor of 0.95 was used to account for a small difference 
in signal level in the two arms. We emphasize that in- 
troduction of this factor lowered the shot noise reference 
level and only reduced the magnitude of the non-classical 
effects we describe below. 

The results of the measurements for Ppump = 
23 mW and 3.3 mW of harmonic power on each out- 
put are shown in Fig. |H1 At 5 MHz, the sum vari- 
ance was (0.50±0.15) dB below the shot noise level 
and the normalized noise variance (amplitude corre- 
lation) was (A|i+ 1)2=1. 78±0. 07. The difference vari- 
ance was (0.10±0.15) dB below the shot noise level giv- 
ing a normalized noise variance (phase correlation) of 
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FIG. 6: (color online) Noise spectra measured by detectors 
Da, Db corrected for electronic noise. The intensity sum which 
gives the amplitude quadrature (red curve) and intensity dif- 
ference which gives the phase quadrature (blue curve) are 
shown with the shot noise level (A|ii2|)^ (black curve). The 
inset shows the sum and difference spectra relative to the mea- 
sured shot noise level together with the corresponding theory 
curves calculated with the same parameters as in Fig. |S|but 
with Ppump = 23 mW. The spectra in Figs. I6I7I were calcu- 
lated with a numerical resolution bandwidth of 200 kHz. 




4- 5 6 r 

Frequency (MHz} 



FIG. 7: (color online) Noise power of the sum of noise vari- 
ances (red curve) together with the quantum noise limit 
2(A|ii2|)^ (black curve). The inset shows the difference of 
the two curves together with the theoretical result using the 
same parameters as in Fig. |S| 



(A|i_ 1)2=1. 95±0. 07. We find 

{A\i+\f + (A|i_|)2 = (1.78 ± 0.07) + (1.95 ± 0.07) 
= 3.73 ±0.14 < 4, 

which satisfies the inseparability criterion given in Eq. 
(0). Expressed in dB the entanglement signature is 
101ogio(3.73/4) = -0.3 dB. 

The sum of the measured noise variances (Ajz+D^ -|- 
(A|i_|)2 is shown in Fig. together with the quantum 
noise limit (QNL) of the sum of the amplitude and phase 
quadratures of the beams which is given by 2(A|zi2|)2 for 
coherent state beams. We see that the sum of the noise 
powers is below the QNL which represents an experimen- 
tal demonstration of non-separability of the two beams. 
The degree of inseparability is close to the theoretical 
prediction at low frequencies and falls at high frequency 
which we attribute to a roll-off in the detector response. 

In conclusion we have demonstrated a new approach 
to generating bright entangled beams using SHG in a 
cavity with two output ports. The amount of entangle- 
ment observed was limited by technical issues, particu- 
larly the available pump power. Calculations [llj show 
that (A|i+1)2 + (Ali_l)2 ~ 2.4 should be possible using 
available nonlinear materials with Ppump ^ 0.5 W. This 
technique may be of particular interest for applications 
such as quantum lithography 3J that will benefit from 
short wavelength entanglement, and could also be a con- 
venient source of tunable entanglement at atomic transi- 
tion frequencies in the visible and near infrared part of 
the spectrum if a longer wavelength infrared pump laser 
is used. 
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